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Abstract 



In this paper we analyse the exact solutions to scalar PDEs obtained 
thanks to summable Taylor series provided by Adomian's decomposition 
method. We propose a modification of the method which makes the cal- 
\ culations of Taylor coefficients easier and more direct. The difference is 

' essential for instance in case of nonhomogenous equations or initial condi- 

tons and is illustrated by some examples. 
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1 Introduction 

The area of exact solutions to nonlinear differential equations has become very 
popular in recent decades, when the development of personal computers enabled 
more efficient work with known algorithms. Adomian decomposition method 
[U El |3] in the matter of fact was developed to find approximated solutions to 
differential equations, but in many publications j3l 2| we can find interesting ex- 
amples where obtained power series were actually summable to exact solutions. 
Typical way to obtain such solutions is to sum up certain Taylor series. In our 
paper we are going to present some situations when it seems reasonable to use 
modified techniques to obtain Taylor series. In mathematical physics we often 
have to deal with scalar PDEs of space and time variables x, t 6 R. We will show 
that for nonautonomous equations of this type the method could be easily mod- 
ified to get Taylor series directly. In order to explain our idea let us first briefly 
present the classical method, so it would be easier to show the differences. 
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2 Adomian decomposition method 

Let us consider the following one- dimensional equation in Cauchy-Kovalevska 
form: 

u t (x, t) = F(u, u x , u xx , ...,u x n) + g(x), (1) 

where x,teR,u t = ff, u x , = 0, F(p ,p 1 ,p 2 ,P3, -,p n ) is a polynomial function 
of it's arguments and g(x) is an analytic non-autonomous term. The correspond- 
ing initial condition is u(x, 0) = fo(x). This quite special case would be perfectly 
sufficient to present the advantages of proposed modification. 
Let us introduce the auxiliary notation: G[u] = G(u, u x , u xx , u x n) for any func- 
tional G defined on jet-space. In classical approach LHS of ([T|) is usually split into 
two parts: F[u] = Lp[u] + N F [u], where L F [u] is a linear operator with respect 
to while A^^[w] is nonlinear part of F[u}. Then the operator 

L~\.)= f\.)dt 
Jo 

can be introduced to express the solution of (JTJ in the form: 

u = fo( x ) + g( x ) t + / L F [u] + N F [u] dt. 
Jo 

Next we assume u = Yli^o Ui anc ^ consequently Lp[u] = Yli^o Lf[ui] and Np[u] = 
^'Ei^o] = Si^o^*- ^ ne new ly introduced terms A* are so-called Adomian's 
polynomials, which could be obtained e.g. with the help of following formulae: 

n i=0 

In the paper [6] author present very intuitive way to obtain these polynomials. 
The idea could be easily understood from the example below. 

Example 1. For instance if we take nonlinearity in the form Nf[u] = uu x then 

N F [u] = (u + eui + e 2 u 2 + ...)(u Qx + eu lx + e 2 u 2x + ...) = 

= Uo u 0x + + u u lx ) + ... 

and here Ai would be a coefficient at e\ 

Let us notice that "e"-notation was not used in original paper, but in our 
opinion it makes the choice of Ai more clear. Let us underline that each polynomial 
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Ai is dependent only on the functions u ,...,Ui. No higher orders are involved. 
Going back to the decomposition algorithm: 



u 



f (x)+g(x)t + / Lf[ui] + Ai dt, 

Jo ,_„ ,_„ 



i=0 i=0 

therefore the following recurrence could be defined: 



u 



uo = fo(x) +g(x)t 

ui— L f [uq] + A dt 
Jo 

= / L F [u n -i] + A n _idt 
Jo 



k 

U 



lim N Ui. 



fc— >00 ' 

i=0 

Remark 1. If g(x) = i/ien there exists a sequence of functions hi(x) such as 
u n = h n (x) t n . 

Proof. By induction uq = fo(x) =: ho(x). Let us assume, that Ui = hi(x)tf so t 
appears in the same power as e in the example [1] which means that A, = ki(x) t l 
and 

u i+1 = / L F [hi(x) f ] + kiix^dt = - — -(L F [hi(x)} + h(x)) =: h l+1 t i+1 . 

i + 1 



o 



□ 



Remark 2. If fo(x) = £/ien t/iere exists a sequence of functions hi(x) such as 

U n = h n {x) t n+1 . 

Proof. By induction u = t g(x) =: h (x). Let us assume, that itj = hi(x)t l+1 
which means that Ai = ki(x) t l+1 and 

u i+l = [ L F [hi(x) t i+1 ] + h(x)t i+1 dt = ^-(L F [hi(x)] + ki(x)) =: h l+l t l + 2 . 



□ 

These remarks imply that in case of autonomous equation or zero initial con- 
dition the algorithm leads straight to Taylor series. 
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3 Main Results 

In our following research it would be comfortable to skip dividing F[u] into two 
parts. The whole functional F[u] could be as well approximated by Adomian 
polynomials. Let us also denote by = Yl^o^- ^ ne key difference would be 
the fact that now we choose polynomials in different way. Let B = F[uq], but 

B< = F[u {l) ] - Flu^-V] 

for i = 1, 2, ... . Comparing to the example [1] now we obtain B = A = Uq Uq x , 
but for a change B\ = F[w + u±\ — F[uo] = Uo U\ x + U\ Uq x + U\ u± x . However it 
is obvious that still lim^oo Bi = F[u]. The new recurrence is: 



uo = fo{x) +g(x)t 

ft ft 

F[u }dt= / B dt 
Jo 



u 2 



[ B 1 dt= [ F[u Q + Mi] - F[uq] dt 
Jo Jo 



Uk 



[ F[u {k -^]-F[u^]dt. 
Jo 



Theorem 1. Let us consider a partial differential equation in the form u t (x, t) = 
F[u(x, t)]+g(x), together with the initial condition u(x, 0) = fo(x) where x,t el, 
u : R 2 — > R, F[u(x,t)] = F[u(x, t), u x (x, t), u x n(x, t)], u x i = We also 

assume, that F[u(x,t)} is analytical of it's arguments and F[0] = 0. Then formal 
Taylor series for the solution u(x,t) could be found using formula (TJ|). 

Proof. Let us start with the summation: 



u 



(k) 



Is . 

J2ui = u + [ F[u^]dt. 
i=o Jo 



Therefore using Taylor's formula 

< (fc+1 > = u Q + [ F[vF>] dt 
Jo 



u + f F[u^]\ t=Q - 
Jo 



dF[ U M] duW + dF[u^} du x k) 







9( M W) 2 


\ dt ) 



duW dt du W 

8 2 F[u^] du^ du ( x k) 
+ 2 du( k )du ( x k) ~dt dF 



dt 



... + 



dF[vW\ du x k } 



dt 



t+ 



J t=0 



dF[uW] d 2 u (k ) 



t=o 



-+.... dt 
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u 



+ F[u<-% =Q t + 



dF[u^} du^ dF[uW] du { x k) dF[u^} du {k } 



~d 2 F[ U W] 




0(u(*O) 2 


{ dt J 



duW dt du ( x k ) dt 
d 2 F[u^} du^ du x k) 



+ 2 



du^du x k) dt dt 



" du ik ) dt 
dF[u^\ d 2 u (k ) 



t=0 
.3 



2! 



t=0 



t' 

3! 



Before we continue let us formulate the following lemma: 

Lemma 1. If uW = a + ait + a 2 t 2 + ... and = b + b\ t + b 2 t 2 + ... then 

a s = b s for s < i. 

Proof. The statement holds if and only if Ui + \ = t l+1 hi+i(x, t) for some analytical 
function h k and k > 1. 

Basis: for k — 1 U\ — J * F[u ] fit. Since F is analytical, it could be written in 
series form: 

F(u,U x1 ... } U x n^ ^ ^ t>ioii •••'^x n ' 

io,h,...,i n 

Thus F also could be written as series with respect to t, F[uq] = Cq + C\ t + ... 
and after integration u\ = c 1 + % t 2 + ... . 
Inductive step: we assume, that u k = t k h k (x,t), then 

u k+1 = [ F[u +u 1 +...+u k }-F[uo+ui+...+u k - 1 ] dt = [ F[S+t k h k {x,t)]-F[S]dt, 
Jo Jo 

where S = u + U\ + ... + u k -\. Using the series form: 
u k +i = J b ioil .. An ((S+t k h k (x,t)y°(S+t k h k (x,t))i\..(S+^ 

In the main theorem we assumed F[0] = so &00...0 = 0. Thus the smallest 
possible power of t in the sum is t k and after integration the proof is completed. 

□ 



Going back to the main proof: 
u Q +F[fo] t+ 



f[/o](9W + F[f ]) + §^[fo](9(x) + F[f ]) x + ... + ^-{f Q }(g( x ) + F[f ]) x , 



0[/o] (9(x) + F[f ]) 2 + ... + 2^-[f ](g{x) + F[f ])(g(x) + F[f ]) x + 



dF r ^fdF r ^,,^ _ r „ 1N dF,^, /s „ r „„ dF \ 1 t 



+ g^iM [faVoM*) + F IM) + q^VoMx) + F[f ]) x + ... + —[f ](g( X ) + F[f ]) x 
In other words, after taking limit we obtain the following formal form: 

4-2 



t t 
u{x, t) « a + ai — + a 2 — + ... 



(2) 
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where 

a o = fo 

oi = g(x) + F[oq] 

dF dF dF 

a-2 = -Q^[ a o] a i + g^~\- a °\ a ix + ■■■ + -q^— [ao] a>ix n 

d 2 F d 2 F d 2 F dF dF 

«3 = ^Na?+^MaL+...+2^Wa 1 a 1 , + ... + -[/ ]a 2+ ...+ — [/ ]a 2a 



Further terms could be easily obtained using formulae for higher differentials and 
the key fact, that u k and u h+1 have the same coefficients up to kth power. □ 



4 Examples 

Example 2. Let us start with the example from ^j. Authors considered the fol- 
lowing Fisher's equation 

u t = u xx + 6u(l-u) (3) 
together with the initial condition 

1 

u(x, 0) 



and obtain the exact solution of |3|) using Adomian decomposition method. Here 
F[u] = F(u,u xx ), d -^l = 6 - 12 «, 9 -^^ = I, = 1, /„ = jj^, 

a o — fo 



(l + e x ) 2 J (l + e*) 2 ^ (1 + e x ) 2 J (1 + e x f 

e x (2e x - 1) 



a 2 = (6 - 12/ )ai + ai xx = 50- 



a 3 = -12(ai) 2 + (6 - 12/ )a 2 + a 2ra = 250 



e x + l) 4 

4e 2x _ 7e x + x 



v e x + l) 5 
etc.. 



The result after summation repeats the cited paper: 

- 1 
u[x,t)- ^ + eX _ m)2 - 
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As we have seen from propositions [TJ [2] the classical method does not lead 
straight to Taylor's series only if the equation or initial condition is nonhomoge- 
nous. The next two examples cover both situations. 

Example 3. Now let us consider the non- autonomous heat equation f^: 

tk = u xx + sinz, 

with the initial condition u(0,x) = cosx. Using classical approach we obtain: 

u = cos x + t sin x 
1 2 • 

U\ = — tcOSX -t SIM 

2! 

1 2 Is- 
U2 = —jt COS X + —t Sill X 

Applying theorem[7\we can directly obtain Taylor series (the only non-zero deriva- 
tive is -J=^— ): 

a = cosx 

ai = sin x + [cos x] xx = sin x — cos x 
dF 

a 2 = — a\ xx = — sin x + cos x 

(JU XX 

dF 

a 3 = — a 2xx = sin x — cos x 

ou xx 

etc., 

so finally 

t 2 t 3 

u(x,t) ?s cosx + (sin x — cosx) (t H — - — ...) = cos xe * + sinx(l — e ). 

2 . 3 . 

To complete the illustration we choose the example with nonlinearity and 
nonautonomous term. 

Example 4. The following inhomogeneous advection problem is solved in f3j/: 

u t + uu x = x, u(x, 0) = 2. 

With the help of decomposition method the following recursive relations was ob- 
tained: 

Uq = 2 + x t 
ui = -t 2 - ^xt 3 

■ * 
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u 2 = —t 4 H 

2 12 15 

Meanwhile using theorem^ 

fa = 2, g(x) = x, F u [u] 

F Ux [u] = -u F Ux [2] 
and all other derivatives vanish in f . 

a = 2 

a\ — x 
a 2 = -2 
a 3 = — 2x 
etc.. 

In both methods we obtain 

2 ( 1 - h 2 + ^ t*+ ..) +x(t--t 3 + -A = 
V 2! 4! J V 3 15 / 

= 2 sech t + x tanh t. 
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